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Abstract. 

Multi-particle non-equilibrium dynamics in two-channel asymmetric exclusion 
processes with narrow entrances is investigated theoretically. Particles move on two 
parallel lattices in opposite directions without changing them, while the channels are 
coupled only at the boundaries. A particle cannot enter the corresponding lane if the 
exit site of the other lane is occupied. Stationary phase diagrams, particle currents 
and densities are calculated in a mean-field approximation. It is shown that there 
are four stationary phases in the system, with two of them exhibiting spontaneous 
symmetry breaking phenomena. Extensive Monte Carlo computer simulations confirm 
qualitatively our predictions, although the phase boundaries and stationary properties 
deviate from the mean-field results. Computer simulations indicate that several 
dynamic and phase properties of the system have a strong size dependency, and 
one of the stationary phases predicted by the mean-field theory disappears in the 
thermodynamic limit. 
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1. Introduction 

Asymmetric simple exclusion processes (ASEP) have become an important tool for 
understanding complex non-equilibrium phenomena in chemistry, physics and biology 
[H [2] . In the absence of fundamental framework for non-equilibrium systems ASEP have 
emerged as reference models, and they play the role similar to the Ising model in the 
studies of equilibrium systems [3J . ASEPs have been applied successfully to analyze the 
kinetics of biopolymerization [1], protein synthesis [SlEllTj, molecular transport through 
nanopores and channels [8], the motion of motor proteins along cytoskeleton filaments 
P, [ini [H], and investigation of car traffic processes [I2l [13]. 

ASEPs are one-dimensional models where particles, that interact via an exclusion 
potential, hop along discrete lattices. Many non-equilibrium processes can be described 
by single-channel exclusion processes. However, the necessity to analyze more realistic 
complex phenomena, that involve the transport along the parallel channels and existence 
of internal states, e.g., motor proteins, vehicular traffic and hopping of quantum 
dots, stimulated the development of multi-channel ASEPs [HI [T5| [T6l [T7t [T8l [T9| [20] . 
Theoretical analysis of multi-channel exclusion processes indicates that a coupling 
between channels has a strong effect on stationary-state phase diagrams and it leads to 
several unusual phenomena, such as localized domain walls [T9l [20]. 

One of the most intriguing phenomena observed in ASEPs is symmetry breaking 
when the microscopic symmetric dynamic rules lead to the existence of macroscopic 
asymmetric stationary-state properties for some sets of parameters [TH [2T1 [22l [231 [Ml [251 
[26], [27|, [28] . However, mechanisms of symmetry breaking phenomena in ASEPs are still 
not well understood [25]. Currently, these phenomena are mostly described by mean- 
field approaches that neglect the correlations in the system, however, the agreement with 
computer simulations results is rather qualitative [211 [221 [231 [2l] . There are suggestions 
that dynamic amplifications of fiuctuations might lead to a broken symmetry in the 
asymmetric exclusion processes [28] . It is also argued that the origin of this phenomena 
might be due to the effective boundary defects [25j. Note also that symmetry-breaking 
has not been obtained in the exact solutions of ASEPs [221 129] . 

Originally, the symmetry breaking has been observed in the single-lane exclusion 
process with two particles moving in opposite directions [211 [22|- Using a mean-field 
approach and Monte Carlo computer simulations, it was shown that two stationary 
phases with broken symmetry are possible for some part of the parameter space, 
although, surprisingly, the agreement between theoretical predictions and numerical 
calculations was not good, in contrast to other steady-state properties of other ASEPs 
where the mean-field description is rather successful [H [2]. There are controversial 
reports about the number of phases with broken symmetry. The existence of one of 
the symmetry-broken phase has been disputed on the basis of computer simulations by 
Arndt et al. p3]. More extensive computer simulations [24J confirmed the presence of 
two asymmetric phases, but also it was argued that phase transitions to the symmetric 
phase is different from the mean-field predictions. However, extensive high-precision 
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Monte Carlo data pZj indicated that the disputed intermediate phase is a finite-size 
effect and it disappears for very large systems. It is not clear if this observation of 
single asymmetric phase is correct for all ASEPs where the symmetry breaking have 
been observed. In addition, in single- channel two-species ASEP, where the symmetry 
breaking is observed, particles of different types interact with each other at every site 
of the lattice, i.e., the interaction is global. It is not clear how the interaction between 
these particles localized in the specific parts of the system will affect the symmetry 
breaking. 

In this paper we investigate two-lane ASEPs with two types of particles moving 
along different channels in opposite directions. The particles cannot jump between 
the lanes, and two channels interact with each other only at the boundaries, the 
problem is motivated by the cellular transport of kinesin and dynein motor proteins 
that move along microtubules in opposite directions [30j . The system is analyzed using 
a mean-field theoretical approach and extensive Monte Carlo computer simulations at 
different system sizes, and we conclude that the symmetry breaking is taking place in 
this system for some range of parameters. This observation supports the idea that 
symmetry breaking can be viewed as an effective boundary-induced phenomenon. 

The paper is organized as follows. In section 2 we give a detailed description of 
the model and we solve it in the mean-field approximation. In section 3 we present 
and discuss Monte Carlo computer simulations, and we compare them with theoretical 
predictions. Finally, we summarize and conclude in section 4. 

2. The Model and Theoretical Description 
2.1. Model 

We consider a system of two parallel one-dimensional lattices with two types of particles 
moving in the different lanes in the opposite directions, as shown in Fig. 1. Two 
channels are identical and both have L sites. The particles are interacting with hard- 
core exclusion potential which means that each lattice site can be empty, or it can be 
occupied by no more than one particle. The hoppings between the channels are not 
allowed, and the particle in the channel 1 (2) moves to the right (left) with the rate 1 
if the next site is empty. The particles can enter the corresponding channel with the 
rate < a < 1 if the first site at the given lane is empty and the site L at the other 
channel is unoccupied: see Fig. 1. The particles exit both lanes with the same rate 
< /3 < 1, that does not depend on the occupancy of the site 1 in the other channel. 
Note that the particle at site L — 1 can proceed to the site L if the exit site is unoccupied 
independently of the status of site 1 at the other lattice. Thus this model can be viewed 
as two-channel ASEP with narrow entrances. 

The model, as described above and illustrated in Fig. 1, is similar to a single- 
channel ASEP with two species where symmetry breaking has been previously observed 
[2T| [22| [23| [211 [27j. However, there are several important differences between two 
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Figure 1. Schematic view of two-lane ASEP with narrow entrances. Dark-gray 
particles move only along the channel 1, while light-gray particles can be found only 
on the channel 2. Allowed transitions are shown by arrows. Inter-channel transitions 
are not allowed. Entrance rates at both channels are equal to a if there are no particle 
at the exit of the other lane. Exit rates are equal to /3 at both channels independently 
of the occupation status of the exit site at the other channel. 



systems: in our model there are no restrictions on the bulk densities of particles in 
both channels and only the entrance to the channel is influenced by the other channel. 
The interaction between two types of particles is local in our system, in contrast to 
single- channel two-species ASEPs where particles interact along the whole lattice. 



2.2. Mean-Field Theory 

Theoretical studies of ASEPs indicate that an approximate mean-field analysis provides 
a very good description of these processes, as confirmed by comparison with available 
exact solutions and computer simulations [H [2] . Following this approach, we introduce 
Pi and nii as occupation variables for the channel 1 and 2, respectively, such that Pi = I 
and/or = 1 if the corresponding site i in the channel 1 and/or channel 2 is occupied, 
and otherwise. The mean- field theory then assumes that there are no correlations in 
the probability of finding two particles at any two sites, i.e., 

{PiPj) = {Pi){Pj) = PiPj, {mirrij) = {mi){mj) = uiimj. (1) 

Our system with particles moving on two channels in opposite directions can be 
viewed as two single-lane ASEPs with exit rates 13 and entrance rates ai and 0^2 for the 
lane 1 and 2, respectively. From the dynamic rules of the system we conclude that the 
effective entrance rates are given by 

ai = a;(l - mi), 0:2 = «(! - Pl)- (2) 

Large-time dynamics and properties of single-channel ASEP are known exactly [HE]. 
In this case there are three stationary phases, specified by the processes at the entrance, 
at the exit and in the bulk of the system. For a < (3 and a < 1/2 the entrance is the 
rate-limiting step in the overall dynamics, and the system is found in a low-density (LD) 
phase with the particle current and bulk density 

JlD = Ol{^- Oi), Pbulk,LD = Ol- (3) 
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When the exit is the slowest process in the particle dynamics (for a > /S and (3 < 1/2), 
the system is in a high-density (HD) phase with steady-state properties given by 

JhD = P{1-P), Pbulk,HD = I - P- (4) 

For a > 1/2 and (3 > 1/2 the dynamics of the system is determined by processes in the 
bulk, and we have a maximal- current (MC) phase with 

JmC = ^5 Pbulk,MC (^) 

Since every channel can be found in one of three stationary phases, there are 
nine possible phases in the two-lane ASEP with narrow entrances, assuming that the 
dynamics in the channels is independent of each other. However, some of the phases 
might not exist. Because of the coupling at the boundaries of the channels some of the 
phases could be asymmetric. To investigate this possibihty let us consider separately 
symmetric and asymmetric cases. 

2.2.1. Symmetric Phases In symmetric phases all properties in both channels are 
identical, 

Jl = J2, Pi = P2, "1 = "2, (6) 

where J, and pi are the stationary-state particle current and the bulk density in the 
channel i (with i = 0, 1). Because of the symmetry, the probabilities to find particles 
in both channels are related via 

Pi = ITLL-i+l. (7) 

In MC phase the stationary current is equal to Ji = J2 = 1/4. From the expression 
for the current at the exit of the channel 2, 

J2 = /3mi = 1/4, (8) 

we derive mi = Substituting this result into Eq. ([2]) yields the value for the effective 
entrance rate, 

«i = «(1 - (9) 

Then the conditions for the existence of this phase (ai > 1/2 and (3 > 1/2) can be 
written as 

In LD phase the current in both channels is equal to Ji = J2 = ai(l — Oii). The 
exit current from the lane 2 is equal to 

J2 = /3mi = ai(l - ai), (11) 

and we obtain mi = _ After substituting this result into Eq. ([2]), it can be shown 

that 



a + (3- J(a + py - 4q;2/5 

= —a • (^2' 



Symmetry Breaking in Two- Channel Exclusion Processes 



6 



Because in the low-density phase ai < P and ai < 1/2, one can show that the condition 
for existence of this phase is the following. 

Another possible symmetric phase is HD phase where Ji = J2 = f3{l — (3). Because 
the exit current in the channel 2 is J2 = f3mi, we can show that mi = 1 — j3. The 
existence of this phase is specified by the condition ai = a(l — mi) > (3, and it can 
take place only for a > 1. This suggests that the HD phase cannot be found in two- lane 
ASEP with narrow entrances. 



2.2.2. Asymmetric Phases In asymmetric phases the currents and densities in each 
channel are generally different, 

Ji ^ J2, ai ^ a2. (14) 

First, let us consider the possibility of phases when one of the lanes supports the 
maximal-current phase. For convenience, assume that the channel 1 is in the maximal- 
current state. The MC/HD phase is specified by 

f ai>l/2 and /3> 1/2, 

[ a2> P and /3 < 1/2. ^ ^ 

However, these conditions contradict each other, and we conclude that MC/HD phase 
does not exist. 

Similarly, the conditions of existence for MC/LD phase can be written as 

r ai > 1/2 and (3 > 1/2, 
\ a2 < /3 and 0:2 < 1/2. 

The stationary currents in the system are given by 

Ji=/9pL = l/4, J2 = /5mi = a2(l-a2), (17) 
which leads to the following expressions, 

1 «(1 - ^) [1 - (1 - Jp\ 



PL = 



mi 



4/3' p 
Because cti = a{l — mi) > 1/2, it requires that 



a 

J 



2a ^ Ap' 



> 1. 



(18) 



(19) 



However, this expression cannot be satisfied for < a < 1 and < /3 < 1, and it leads 
us to conclusion that the MC/LD phase also cannot exist in the two-lane ASEP with 
narrow entrances. 

The situation is different for the HD/LD phase which is defined by the following 
expressions, 

j ai > P and P < 1/2, 
\ a2 < P and a2 < 1/2. 



(20) 
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The stationary current in this phase can be written as 

Ji = PpL = P{l-P), J2 = /?mi =a2(l-«2), (21) 

and these expressions lead to 

ai = a{l — a + a'^P), a2 = ap. (22) 

Combining these equations with the conditions (120| ) we derive the following condition 
for the existence of this phase, 

P < q (23) 

Another possible asymmetric phase is the LD/LD phase, where both channels are 
in the low-density state but with different particle currents and bulk densities. Formally, 
this phase can be described by 

{ ai < P and ai < 1/2, 
\ a2 < P and 02 < 1/2. 
The stationary currents in both channels are given by 

Ji = PPl = otii)- - Oil), J2 = p-mi = 02(1 - 02) , (25) 
and from these expressions we derive 



(24) 



ai = a 



1 - 



02(1 — ^2) 



P 



1 - 



ai(l — Oi) 



P 



a2 = a 

To simplify calculations, assume that ai > 02 and introduce two auxiliary functions 

S = ai + ^2, and D = ai — 02- 
From Eqs. ( l26l) one can show that 



D = — («! - aj- a2 + a\) = -^-0(1 - S). 

Since in the asymmetric phase D 7^ it leads to 

S=l-p/a. 
Similarly, from Eqs. fl26l) we obtain 



S 



a r 



S-{S^ + D^)/2 



Then, using the result (l29l) . it can be shown that 
D 



\ 



1-4/3 + 2 



1 
a 



a . 



The effective entrance rates can be easily calculated in this case from Eqs. 

(EH), 



ai = {S + D)/2 



(26) 

LS, 

(27) 
(28) 
(29) 
(30) 



(31) 
and 



a2 = {S-D)/2 = - 




(32) 
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This phase can only exist if the expression inside of the square root is positive. This 
requirement, in addition with Eqs. fl2^ . determines the conditions for existence of the 
LD/LD asymmetric phase in the mean-field approximation, 



The resulting phase diagram, obtained via the mean-field calculations, is shown in 
Fig. 2. Similarly to the single- lane two-species ASEP where the symmetry breaking 
was first observed [211 1221 [231 [211 [27], the region for existence of the LD/LD asymmetric 
phase is very small. But there is a MC phase in our system, not observed in the 
single-lane two-species model. In addition, the mean-field approach predicts a first-order 
phase transition between HD/LD and LD/LD asymmetric phases, and continuous phase 
transformations for boundaries between asymmetric LD/LD and symmetric LD phases, 
and between symmetric LD and MC phases. 

3. Monte-Carlo Simulations and Discussions 

Our theoretical approach neglects correlations between particles, and it is necessary to 
test theoretical predictions with computer simulations. We performed extensive Monte 
Carlo computer simulations to analyze two-channel two-species ASEP with narrow 
entrances. Theoretical calculations assume that the system is in the thermodynamic 
limit, i.e., L cxo. However, in computer simulations only the processes with finite 
size channels can be tested. The systems with lattice sizes ranging from L = 100 up 
to L = 12, 000 have been used in our simulations, but most computations have been 
done for L = 1000. To accelerate the simulations we implemented the so-called BKL 
algorithm [31]. This method skips uneventful Monte Carlo steps with simultaneous 
correction in the calculation of the time. Number of effective Monte Carlo steps per 
lattice site in our simulations was typically between 2x10"^ and 5x10^. Usually, first five 
percents of steps have been omitted from the calculation because we are only interested 
in the stationary state of the system. 

Since the phase diagram predicted by the mean-field theoretical approach for two- 
channel ASEP with narrow entrances is very complex, it is not enough to measure in 
the simulations only the average densities and the currents per channel. It was proposed 
earlier [231 [21] to utilize also density distribution functions that measure the frequency 
of simultaneously observing average densities for different species. It was shown that 
this method was powerful enough to map precisely the phase diagram of the single- 
lane two-species ASEP. In our system the density distribution function depends on the 
average densities in the first and the second channels determined at the same time. To 
construct the corresponding distributions these densities have been measured after every 
L/10 Monte Carlo steps. The resulting histograms are shown in Fig. 3 for the fixed 
value of a = 0.9, where two asymmetric and two symmetric phases can be clearly seen. 

In addition, to determine exactly the position of the boundary between symmetric 
LD and MC phases the derivative of the average current as a function of the exit rate 




(33) 




nl ■ 1 1 1 ■ 1 ■ 1 

T).2 0.25 0.3 0.35 0.4 

[b] P 

Figure 2. Phase diagram for a two-channel ASEP with narrow entrances. Lines are 
predictions from the mean-field theory. Thick solid lines represent first-order phase 
transitions, while thin solid lines correspond to continuous transitions. Symbols are 
from Monte Carlo computer simulations for the channels with L = 1000. Filled squares 
describe the continuous transitions between symmetric LD and MC phases. Filled 
diamonds describe the continuous transitions between symmetric LD and asymmetric 
LD/LD phases, while unfilled diamonds are at the position of first-order phase 
transition between HD/LD and LD/LD asymmetric phases, a) Phase diagram for 
all parameter space; b) phase diagram for 0.2 < /3 < 0.4. 



(3 has been computed for different fixed values of tlie entrance rate a. In tlie MC phase 
the current reaches a constant maximum value, and we associate the phase boundary 
with the value of /3 when the derivative reaches zero or starts to fluctuate around zero. 
This procedure is very robust and it is illustrated in Fig. 4. 

One of the most controversial issues in studies of symmetry breaking in the single- 
lane two-species ASEP is the question of existence of the asymmetric LD/LD phase 
[23t EH [27] . Similar situation is also observed in the two-lane two-species ASEP with 



Symmetry Breaking in Two- Channel Exclusion Processes 



10 



<x=0.9, p=0.245 



oc=0.g, p=0.255 



[a] 



[b] 




[c] 




oeO.9, p=0'2595 




Figure 3. 3D plots of the probability density P{pi,p2) as a function of the average 
particle densities in both channels for L = 1000 and a = 0.9. Contour plots are 
shown projected onto p\ — p2 plane, (a) asymmetric HD/LD phase at /3 = 0.23; (b) 
asymmetric HD/LD phase at /3 ~ 0.245; (c) asymmetric HD/LD phase at /? = 0.255; 
(d) asymmetric HD/LD phase at /3 = 0.258; (e) phase coexistence between asymmetric 
HD/LD and LD/LD phases at /3 = 0.2595; (f) asymmetric LD/LD phase at /3 = 0.262; 
(g) phase coexistence between asymmetric LD/LD and symmetric LD phases at 
(3 = 0.2685; (h) symmetric LD phase at /? = 0.28; (i) symmetric MC phase at /3 = 0.95. 
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Figure 4. Computer simulations results for (a) current and (b) derivative of the 

current as a fimction of the exit rate /3 for L = 4000 and a = 0.9. Arrows show the 
current saturation point that we associate with the phase boundary between symmetric 
LD and MC phases. 
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narrow entrances. Our mean-field theory predicts that the asymmetric LD/LD phase is 
real, ahhough the range of existence is quite small - see Fig. 2. Monte Carlo simulations 
performed for the channels of size L = 1000 also suggest that this phase is a stable 
stationary-state regime. However, it was shown earlier that there are very strong size- 
scaling effects in ASEPs with symmetry breaking [27]. In order to study this size- 
scaling dependency in our model we carried out computer simulations for the systems 
with different sizes (up to L = 12, 000). Phase boundaries have been computed in Monte 
Carlo simulations for different system sizes by utilizing the density-distribution functions 
and current derivatives method, as described above, and the results are presented in Fig. 
5. It is shown in Fig. 5a that the boundary between the the asymmetric HD/LD and 
LD/LD phases does not depend on the system size, while the region of existence for 
the asymmetric LD/LD phase shrinks constantly with increasing L without reaching a 
saturation. This suggests that in the thermodynamic limit (L oo) the asymmetric 
LD/LD phase does not exists. In addition, the results presented in Fig. 5b indicate 
that the symmetric MC phase also depends on the system size. With increasing L 
the boundary between symmetric LD and MC phases shifts closer to the theoretical 
predictions, but it is not clear if it will reach it in the thermodynamic limit. 



0.28 



0.27 - 



-1 1 1 1 1 1 1 1 1 1 r- 



# HD/LD-LD/LD phase boundary 
▲ LD/LD-LD phase boundary 



LD 



LD/LD 
0.26^9 • 



A 



A 







HD/LD 

_i_ 



2000 4000 6000 8000 10000 12000 
L 



0.9 



0.8 



0.7 



0.6 



0.5, 



[b] 



MC/LD boundary, simulations 
MC/LD boundary, theory 



MC 



LD 



_i_ 



2000 



4000 



6000 



8000 



Figure 5. Phase boundaries as a function of the channel size L. (a) Computer 
simulations results for phase boundaries between the asymmetric HD/LD, the 
asymmetric LD/LD and the symmetric LD phases for a — 0.9; (b) Computer 
simulations results for phase boundary between symmetric MC and LD phases for 
(3—1. Solid line corresponds to the phase boundary calculated in the mean- field 
approximation. 



The phase diagram computed from the Monte Carlo simulations is shown in Fig. 
2, while bulk densities and particle currents are presented in Fig. 6. By comparing the 
results of the numerical simulations with the theoretical predictions we conclude that 
the mean-field approach agrees with simulations only for small values of entrance and 
exit rates. For large values of a and /3 theoretical calculations, although qualitatively 
correct, deviate significantly from Monte Carlo computer simulations. In addition, the 
mean-field approach predicts two asymmetric (HD /LD and LD /LD) and two symmetric 
(LD and MC) phases, however our size-scaling analysis of computer simulations (see Fig. 
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5) shows that asymmetric LD/LD phase does not exist in the thermodynamic hmit. It 
is also interesting to note that two-channel two-species ASEP with narrow entrances 
support the maximal- current phase, in contrast to the single-lane two-species ASEP 
[211 [221 [231 [21 [27]. 

The fact that the mean-field approach cannot fully describe phase diagram and 
stationary-state properties of the two-lane two-species ASEP with narrow entrances 
indicates that correlations are important in the dynamics of this system, similarly to 
the single-lane two-species ASEP [211 iSSl [23l [2ll [27] . However, in our system two 
channels do not interact with each other except at the entrance/exit lattice sites, and it 
is very surprising that this local interaction strongly infiuences the macroscopic particle 
dynamics in the system. This observation supports the idea that symmetry breaking in 
ASEPs might occur only in the systems with the effective boundary defects [25] . 
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Figure 6. Stationary-state properties as a function of /3. Lines are predictions 
from the mean-field theory, symbols are from Monte Carlo computer simulations for 
L = 1000. (a) currents in both channels for a = 0.1; (b) currents in both channels for 
a = 0.8; (c) bulk densities in both channels for a — 0.9. 



4. Summary and Conclusions 

We have investigated two-channel ASEP with two types of particles moving in opposite 
directions along the channels with narrow entrances via the mean-field theory, that 
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neglects correlations in the system, and extensive Monte Carlo computer simulations. 
Our theoretical analysis indicates that there are four stationary-state phases in the 
system, with two of them exhibiting symmetry breaking in dynamic properties. 
However, computer simulations only partially support the mean-field picture. Symmetry 
breaking is found in HD/LD phase, and it is shown that the existence of asymmetric 
LD/LD phase is a size-scaling effect: the phase disappears in the thermodynamic 
limit. In addition, mean-field predictions for particle currents, bulk densities and phase 
boundaries agree with computer simulations quantitatively only for the low entrance 
and exit rates, while in other regions of the parameter space the agreement is mostly 
qualitative. 

In our system particles from the different channels interact only at the entrance 
sites. However, surprisingly, this local interaction still leads to symmetry breaking 
and long-range correlations in the system. We conclude that these phenomena can 
be viewed as a result of putting an effective impurity on the boundary of the lattice. 
It is important to develop a better theoretical description of these processes that can 
explain the mechanisms of symmetry breaking in more detail. It is reasonable to suggest 
that the domain- wall method [32], successful in the determination of mechanisms and 
properties of non-equilibrium processes in ASEPs, is the most promising approach for 
future investigations. 
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